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Motivated by recent developments in magnetic materials, frustrated nanoarrays and cold atomic 
systems, we investigate the behaviour of dipolar spins on the frustrated two-dimensional kagome 
lattice. By combining the Luttinger-Tisza approach, numerical energy minimization, spin-wave 
analysis and parallel tempering Monte-Carlo, we study long-range ordering and finite-temperature 
phase transitions for a Hamiltonian containing both dipolar and nearest-neighbor interactions. For 
both weak and moderate dipolar interactions, the system enters a three-sublattice long-range ordered 
state, with each triangle having vanishing dipole and quadrupole moments; while for dominating 
dipolar interactions we uncover ferrimagnetic three-sublattice order. These are also the ground 
states for XY spins. We discuss excitations of, as well as phase transitions into, these states. We find 
behaviour consistent with Ising criticality for the 120° state, while the ferrimagnetic state appears 
to be associated with drifting exponents. The celebrated flat band of zero-energy excitations of the 
kagome nearest-neighbour Heisenberg model is lifted to finite energies but acquires only minimal 
dispersion as dipolar interactions are added. 

PACS numbers: 75.10.Hk, 75.30.Ds, 75.40.-s, 75.40.Mg 


I. INTRODUCTION 

Long ranged dipolar interactions occur in any lattice 
system of interacting magnetic moments. However, the 
assessment of the relevance of dipolar interactions in 
determining the behavior of magnetic systems has wit¬ 
nessed a recalibration in the recent past. This is largely 
due to the advent of several experimental systems that 
shifted the focus away from purely exchange coupled 
magnets where the dipolar interaction is routinely ne¬ 
glected. 

We can identify at least three broad classes of systems 
which have led to this renewed interest in dipolar in¬ 
teractions. The first are the A 2 B 2 O 7 pyrochlore oxides, 
which most closely resemble conventionally studied mag¬ 
netic system^S. For these, as a result of an interplay of 
crystal field effects, geometry and the specific magnetic 
ions involved, the dipolar interations can be appreciable. 
A second class are nanomagnetic arrayfl, collections of 
nanomagnetic islands arranged in a regular pattern using 
lithography. The magnitude of the moments as well as 
the strength of the dipolar interactions can be tuned to a 
great degree by controlling the dimensions and separation 
of the magnetic islands. These systems are much more 
tunable than the thin film systems studied in the past 
with a view to analysing pattern formation and order¬ 
ing via the dipolar interactiorP. Finally, the last decade 
has seen rapid development of magnetic systems of polar 
molecules and atomic ga ses with large dipole moments 
confined in optical lattice^^. 

Of particular interest is the interplay of dipolar interac¬ 
tions and geometrical frustration. On frustrated lattices, 
an exchange term typically gives rise to a macroscopi- 
cally degenerate yet locally strongly constrained ground 
state manifold, usually lacking conventional magnetic or¬ 


der. This constraint can be thought of as restricting the 
space of states, often in a topologically non-trivial way, 
within which dipolar interactions are to be minimised; 
or, conversely, the dipolar interactions can be thought 
of as lifting the degeneracy, akin to the usual order-by¬ 
disorder physics characteristic of quantum and thermal 
fluctuation^. The combination of exchange and dipoles 
can lead to suprising results, such as in the case of spin 
ic^, where the underlying elementary excitations can 
be seen as doubly gauge chargecP (emergent) magnetic 
monopole^. 

Theoretical efforts to study dipolar spins are several 
decades olcP^i^. An early milestone is the work of Lut- 
tinger and TiszrP^who established that the ground state 
for a simple cubic lattice of dipoles is an antiferromag¬ 
netic arrangement of chains of aligned dipoles. There¬ 
after, Maleev®! found that the long range and anisotropic 
nature of the dipolar interactions can stabilise long range 
order in two dimensional magnets - something that is 
prohibited for short ranged isotropic exchange Hamilto¬ 
nians because of the Mermin-Wagner theorem. Indeed, 
for nanomagnetic arrays and cold atoms in optical lat¬ 
tices the study of two dimensions is particularly rele¬ 
vant. For dipoles on the square lattice the ground state 
likely co nsists of antiferromagnetically aligned ferromag¬ 
netic leg j^^bsi and closely related degenerate stateJ^^. For 
the triangular lattice a ferromagnetic phase has been re¬ 
ported for purely dipolar interactions but it was argued 
that other phases like a 120° phase and striped antifer¬ 
romagnetic phases appear for increasing strength of the 
exchange interactiorP^lEol While there is some agree¬ 
ment about the nature of the low temperature phase for 
several of these systems the precise details of the transi¬ 
tion to those low temperature phases are frequently un¬ 
der debate. The principal reason lies in the subtleties 
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involved in the thermodynamic limit in the presence of 
long ranged (and anisotropic) interactions. 

For some otherwise well-studied lattices, not even the 
dipolar ground state is known. A case in point are 
classical dipolar spins on the kagome lattice, the focus 
of this work. The kagome is perhaps the most-studied 
two-dimensional highly frustrated lattice, for which even 
the low-temperature behaviour of a ’simple’ ne arest- 
neighbour Heisenberg model is remarkably intricatd^il^SI 
We investigate in this paper using a combination of 
Luttinger-Tisza (LT) method, spinwave calculations, nu¬ 
merical energy minimisation and extensive Monte Carlo 
simulations the interplay of exchange and dipolar inter¬ 
actions. We find two distinct low-temperature orderings. 
For weak dipolar interactions we observe 120° three- 
sublattice order with zero net moment; while for strong 
dipolar interactions we find a peculiar ferrimagnetic state 
with continuously varying net moment. Thus we have 
two different three-sublattice kg = (0,0) states at weak 
and strong dipolar interactions (Fig. &)• While our re¬ 
sults for the case of strong dipolar interactions predict 
a finite moment per unit cell as in earlier worlP^, our 
extensive simulations and analytic considerations do not 
support the existence of a disordered non-magnetic sub¬ 
lattice. 

The outline of the paper is as follows. In Section |ll] 
we introduce model and conventions used. In Section 
|III| we present the ground-state phase diagram from the 
Luttinger-Tisza approaclP^. This method fails in the 
case of strong dipolar term and hence in Section |IV| we 
perform a numerical search for the ground state and in 
Section [V] we confirm that this state is locally stable 
via a spinwave analysis. Finally, in Section |VI| we an¬ 
alyze our model using a parallel-tempering Monte-Carlo 
method which confirms our predictions for the ground 
states and provides the details of the corresponding rich 
finite-temperature phase transitions. We close with a 
discussion section. 


II. MODEL 


The kagome lattice given in Fig. [^is an Archimedean 
lattic^^, a triangular lattice of triangles. The positions 
of the triangular Bravais lattice points are denoted by R* 
while each site in the unit-cell is labeled by r^, so that a 
site is labeled by R^- = (R* +^j)- Throughout the paper 
the lattice constant is set to 1/2 such that the full 
translation of the three site unit cell is the unit of length. 

The general Hamiltonian of the system of N spins is 


H=Y. ( 1 ) 


Jtf (R) = 2 K 


|R|3 |R|5 


( 2 ) 


Here R/j is the vector between two interacting classical 
0(3) spins S'f(R*’) and S'^(R^), of unit length, k and 



FIG. 1. Kagome lattice (left) and the ground-state phase di¬ 
agram of the model consisting of nearest-neighbour exchange 
J = sin0 and dipolar interactions of strength D = cos 6 
(right). Ai, A 2 are the basis vectors and a dark green tri¬ 
angle denotes the unit cell of three sites. The system exhibits 
long-range 120° order for 7r/2 > 9 > 9i with 9\ — 10.01° 
and ferrimagnetic order for —7r/2 <9 <9^ where 9^ = 1.03°. 
The latter has two spins inclined with respect to one of the 
unit-cell edges by angle ±(j){9). The area between two phases 
possibly contains an incommensurate intermediate regime. 


I index the unit cell, while i and j run over the sites of 
the basis in the unit cell and Greek a and /3 denote the 
components of the vectors x, y and z. The first term 
of the interaction matrix ([^ is the nearest-neighbor ex¬ 
change while the second is the dipole interaction, with 
the nearest-neighbor distance, included for normal¬ 
ization. A factor ^ has been included to avoid double 
counting. J > 0 is the energy scale of the antiferromag¬ 
netic (AFM) nearest-neighbor exchange. The dipolar en¬ 
ergy scale is 


Mo M 
47r R3 ’ 


( 3 ) 


where y is the magnetic moment of the ions. 

We parametrize the relative strength of J and D via 
an angle 9 (Fig. : 

J = sin0, H = COS0, (4) 

with the unit of energy set to = 1. 

Fourier transformation of the Hamiltonian ([^ yields 

J“^(k)5“(-k),Sf(k) (5) 

k,ij 

Jtf (k) = E exp[-*kR^']. (6) 

kl 


We generate the interaction matrix for the dipolar inter¬ 
actions using Ewald summatiorP^, which we confirm by 
the direct lattice summation possible in two dimensions. 
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FIG. 2. Eigenvalues of the interaction matrix Jj“^(k) along 
lines in the Brillouin zone, a) Spectra for dominant exchange 
case, 6 = 7r/4. b) spectra for D = 1, 6 = Q. The eigenvalue 
lowest in energy is generally at the F point ko. 


III. LUTTINGER-TISZA ANALYSIS 



0 

FIG. 3. Inclination angle ±(?i of two of the spins in the unit 
cell (Fig. as function of 6 (Eq. in the ferrimagnetic 
phase. 


We first determine a ground state using the Luttinger- 
Tisza (LT) met ho where it applies. Decomposing the 
interaction matrix into its Fourier components, and de¬ 
noting by Amin(k) the lowest eigenvalue(s) of the inter¬ 
action matrix, we use the fact that the energy of any spin 
configuration satishes the bound 

H > NXra^ni'ko)■ (7) 

If there exists a spin configuration which can be decom¬ 
posed into a linear combination of only the ’optimal!^ 
LT eigenvectors corresponding to these eigenvalues, it is 
a global ground state. This happens if the “strong con¬ 
straint” of unit length for the spins 

|S.P = 1 (8) 

does not conflict with the optimal eigenvectors, which 
however in general have entries with different amplitudes. 
In the latter case, not unusual for non-Bravais lattices, 
non-optimal modes have to be admixed, and the LT ap¬ 
proach only yields an (often rather useful) guess at possi¬ 
ble ground states, or at least at leading instabilities from 
the high-temperature paramagnet. 


A. Dominant nearest-neighbor exchange. 


For pure nearest-neighbor antiferromagnetic exchange 
6 — tt/2 {D = 0), the lowest branch of the interaction 
matrix is exactly flat (di^ersionless) reflecting the high 
ground state degeneracy^BSni Decreasing d we move to 
nonzero D > 0 which immediately lifts the degeneracy, 
selecting a ground state at wavevector kg = (0,0), as 
shown in Fig. a) for 9 = 7r/4. The optimal eigenvec¬ 
tor satisfies the constraint § and results in a 120° state 
which is doubly degenerate reflecting two possible chi¬ 
ralities. Further increase of D leads to level crossing at 
9i = 10.01°. Hence the 120° state is certainly stable up 
to this point, as we have also confirmed in Monte-Carlo 
simulations. 


B. Dominant dipolar exchange. 


For 0 < 01 LT no longer yields an exact ground stat^^. 
Instead, we enter an intermediate regime where neither 
spin-wave nor Monte-Carlo computations (see Sections 
[V| and VI) allow us reliably to conclude on the ground 
state. This regime persists up to the point 02 = 1.03° 
beyond which the 120° state is no longer even a stable 
local minimum at kg = (0, 0). 

For purely dipolar interactions 9 — 0 the minimal 
eigenvalue eg = —2.487 is doubly degenerate and again 
occurs at kg = (0,0), Fig. |b. The best state we find has 
two of the spins are inclined approximately by (j)~ ± 16° 
with respect to one of the unit-cell edges while the third 
spin remains unchanged (right panel of Fig. [^. This sit¬ 
uation persists, with varying (j){9) until the ferromagnetic 
point 0 = —Ti 12. However, in general no combination of 
the pair of eigenvectors satisfies the strong constraint on 
spin length (|^. To determine the true ground state for 
hard unit length spins, we thus need to allow the admix¬ 
ing of other modes, so that we next turn to numerics. 


IV. NUMERICAL ENERGY MINIMIZATION 
FOR D » J 


. Our Monte-Carlo simulations (Section |VI[ ) do unveil 
an kg = (0,0) ordering at low temperatures, suggesting 
that hard spin constraint may optimally be satisfied by 
admixing higher modes at kg = (0, 0) only. We therefore 
constrain our problem to a single unit cell and perform 
a numerical minimization of the Hamiltonian (§. The 
minimal energy configuration for the single unit cell is 
indeed the state found in full lattice Monte Carlo simula¬ 
tions. The ground state is a ferrimagnetic configuration 
in which the spins take the following angles with one 
of the three edges in the unit cell 

'('1 = 02 = —0, 03 = 0, 
















4 




FIG. 4. Spin-wave spectra for a) 120° state 6 — 7r/4 and 
for b) the ferrimagnetic ground state obtained by energy min¬ 
imization at 0 = 0. In a), the lowest branch remains almost 
perfectly flat, while acquiring a gap oc '/D (inset). 



T/S T/S 


FIG. 5. a) Specific heat scaled by the temperature and b) 
temperature dependence of the reduction of sublattice mag¬ 
netization due to quantum fluctuations calculated using linear 
spin wave theory, which is controlled in the limit of small T. 


with 


temperature, e.g. below a scale set by the gap in the 
excitation spectrunP^. 

We first confirm that the 120-degree and ferrimagnetic 
states at 0 = 7r/4, 0, respectively, are stable to quantum 
fluctuations. While it is known from previous studie^^ 
that for the 120° state the spin-wave excitation spec¬ 
trum has a fully dispersionless (flat) band at zero energy 
as well as two-fold degenerate acoustic mode, the ad¬ 
dition of D leads to a gap in the excitation spectrum 
proportional to \fD at small D. We plot the spin-wave 
spectra for cases 0 = 7r/4 and 0 = 0 in Fig. Q where 
it is clearly seen that in both cases the leading effect of 
dipolar interactions is pushing the zero modes to finite 
frequency, expected on account of the absence of a con¬ 
tinuous symmetry. The dispersion of the lowest branch 
of the 120° state is only weakly affected. This fact can 
manifest itself in finite energy almost k-in dependent res¬ 
onance in inelastic neutron scatterin^^nEH existence 
of the gap in the spectrum affects the corresponding spe¬ 
cific heat and sublattice magnetization (Fig. [^. The 
gap leads to an exponential suppression of specific heat 
Cy ~ exp[—A/T] or reduction of staggered magnetiza¬ 
tion AM{T) ^ exp[—A/T] with the gap A. 

Moreover we have checked for both 0 = 0,7r/4 that 
these are the only stable spin-configurations at kp = 
(0,0). We close this section by noting that both the 
120° and ferrimagnetic states are locally stable within 
the boundaries of intermediate phase (Fig. [^. 


(j) « 36.42° 


(9) VI. MONTE-CARLO SIMULATIONS 


As we change —7r/2 < 0 < 02 we can obtain a minimal 
energy ferrimagnetic configuration with a drifting (/(0), 
as shown in Fig. 


V. LINEAR SPIN-WAVE THEORY 


We next study the role of quantum fluctuations around 
the two ground states discussed above. We find that both 
states are locally stable, and exhibit a lowest band with 
little dispersion, in particular for the 120° state. 

We evaluate the spin-wave spectrum of non-collinear 
spin structures using standard method^^^HaH Hamil¬ 
tonian of a Bose gas of magnons reads 


i7 = i7(°)+^^e,(k) (10) 

k i 




a|'(k)aj(k) -|- a|(-k)aj(-k) 


where ai(k) are boson annihilation operators with 
the classical ground state energy. For a stable ground 
state spin configuration, H is Hermitian and all the spin- 
wave eigenenergies ei(k) are real. This yields the spe¬ 
cific heat Cy{T) and magnetization M{T), allowing in 
principle for comparison with experimental data at low 


This section pursues two goals. Firstly, the ground 
states are confirmed numerically; secondly, the corre¬ 
sponding finite-temperature phase transitions are anal¬ 
ysed in detail. This is done with computationally inten¬ 
sive but tractable Monte-Carlo simulations of the system 
on finite lattices with linear dimension of T < 24 unit 
cells or N < 1728 sites. 

We employ parallel tempering with 64 to 128 replicas 
in the temperature range T = 0.125 — 2.95 for the phase 
transition analysis and in the range T = 0.00625 — 2.95 
to investigate the low energy configuration of the dipoles. 
One Monte-Carlo step corresponds to a sweep over the 
lattice in which on average every spin is touched. We 
perform « 10° Monte Carlo steps for the thermalization, 
followed by « 10° steps for every measuring round. 

We obtain thermodynamic properties of the model 
(specific heat, uniform and staggered - 120° state 
-magnetization, magnetic susceptibility, fourth order 
Binder cumulant) as well as the structure of the low- 
temperature spin configuration. 

For the set of parameters leading to the ferrimagnetic 
ground state we analyze the phase transition via the be¬ 
havior of the magnetic order parameter 

^ = ( 11 ) 
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FIG. 6. Monte-Carlo results for specific heat C/N, Binder 
cumulant Ul, magnetization M and susceptibility x for 0 = 0 
{D = 1 and J = 0). 


FIG. 8. Scaling collapse for the magnetic (chiral) order pa¬ 
rameter and its susceptibility for the transition into the 120° 
state (left) and the ferrimagnetic state (right). 



FIG. 7. Monte-Carlo results for specific heat C/N, Binder 
cumulant Ul, chiral order parameter and corresponding 
susceptibility x for 9 = 35.6°. 


where the sum is taken over all the sites in the lattice. 
For the planar 120° ground state order we investigate the 
order parameter which captures the particular chiral spin 
pattern, 


( 12 ) 

where 

H S(R-f rj )exp(f<()^), (13) 

R+Fj 

and (pj are sublattice phase angles 4 >i = 0 , ^2 = 27r/3, 
and (/>3 = 47 r/ 3 . 

To investigate the corresponding finite-temperature 
phase transition we also compute the fourth order Binder 
cumulant 

= ( 14 ) 


as well as susceptibility 

((O^)-(O)"), (15) 

and specific heat per spin 

C/N=~^{{E^)-{Ef). (16) 

To characterize the phase transitions we employ stan¬ 
dard finite-size scaling 

C{L^I'^t) = (17) 

where t = {T — Tc)/Tc is the reduced temperature. To 
obtain the critical exponent l/v and critical point Tc we 
use the scaling relation for the Binder cumulant 

f7(Li/"t) = Ul. (18) 

We extract v, /3, a, 7 and via data collapse. 

Let us first look at the finite temperature transition to 
the 120° state. We perform MC simulations deep in the 
ordered phase for 9 = 35.6°. Collapsing the curves for 
chiral order parameter, Binder cumulant, susceptibility 
and specific heat yields critical temperature Tc as well as 
the full set of critical exponents v, j3 ,x, a see Fig. 
The transition occurs at = 0.692(5), consistent with 
the 2D Ising universality class with critical exponents 
y = I, (3 = 1/8, 7 = 7/4, a = 0, reflecting the discrete 
Z 2 symmetry of the chiral order parameter. 

For dominant dipolar interactions we analyze two 
points, 9 = 0, and for 0 = 1°. The ferrimagnetic or¬ 
der has a six-fold discrete symmetry and at T = 0 differs 
only in the angle (j) of the two inclined spins (Fig. 
b)). We therefore expect the corresponding transitions 
to belong to the same universality class. 
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TABLE I. Critical exponents for the continuous phase tran¬ 
sition analyzed with classical Monte-Carlo. 



l/v 

ajv 

P/v 

7/iz 

U. Class 

6 « 7r/5 

1 

0 

1/8 

7/4 

Ising 

0 = 0 

1.05(3) 

0.10(3) 

0.25(2) 

1.5(2) 

Unknown 

0 = 1° 

1.17(3) 

0.32(3) 

0.25(3) 

1.5(3) 

Unknown 


Our Monte-Carlo data show a clear divergence of the 
ferromagnetic order parameter, specific heat and suscep¬ 
tibility as well as crossings of fourth-order Binder cu- 
mulant curves. This suggests a single second-order phase 
transition from a high-temperature paramagnet to a low- 
temperature ferrimagnetic phase (Fig. [^. Both for 
6 » = 0 and 0 = 1 ° we can extract critical temperatures 
Tc = 0.439(2) and Tc = 0.406(5) as well as the set of 
exponents which lead to the best data collapse of Binder 
cumulant, magnetization, susceptibility and specific heat 
(Table |l|. Note that correlation length exponent u and 
order parameter exponent /? increase monotonically with 
J while the ratios fijv « 0.25 and ^jv « 1.5 remain 
constant with a two-dimensional scaling law implying 
77 = 2/3/7Z « 0.5. 

This appears to provide an example of the so-called 
’’weak universality” hypothesis which states that ratios 
of exponents should be independent of the details of 
system Hamiltonian with 77 = 2(3jv and ^jv univer¬ 
sal while a and /3 are allowed to chang^^. The ’’weak 
universality” behavior is often observed as a drift from 
Berezinskii-Kosterlitz-Thoules (BKT) expone nts to dis¬ 
crete (i.e. Ising, Potts) transition exponent^^HSH^ 
may of course be related to the existence of a large length- 
scale. Note that in our case, we have a correspondence to 
a six-state clock model arising from a Hamiltonian with 
both nearest-neighbor and long-range interactions. Indi¬ 
vidually, a six state clock model with only the former ex¬ 
hibits two KT transitions (not observed here)p2H131 while 
mean-field studies for the case of long range dipolar in¬ 
teractions sugg est a s ingle second order low-temperature 
phase transitioiP^l^. 

A low-temperature phase transition of pure dipoles on 
the kagome lattice was recently observed in the 0{N) 
Monte-Carlo studies in the Ref. US) The nature of the 
low-temperature spin arrangement was however not re¬ 
solved due to the high computational cost of the 0{N) 
Monte Carlo algorithm inversely proportional to temper¬ 
ature. We have investigated the system at significantly 
lower temperatures, where snapshots of the spin con¬ 
figurations give clear evidence of ferrimagnetic order at 
ko = (0,0). At the same time, the temperature depen¬ 
dence of the static structure factor does not indicate any 
intermediate ordering between the low-temperature ferri¬ 
magnetic state and the high-temperature disordered con¬ 
figuration. Together with our LT and spin-wave studies 
this rather strongly suggests that ferrimagnetic k = ( 0 , 0 ) 
state is the low-temperature configuration of the dipoles. 


In the intermediate regime, due to existence of many 
metastable energy minima, our Monte-Carlo simulations 
do not equilibrate even for our extensive parallel setup. 
We thus cannot provide a clear picture of physical quan¬ 
tities and leave a detailed investigation of this possibly 
incommensurate regime for future studies. 


VII. DISCUSSION AND CONCLUSION 

We have determined ground states, excitations, and 
phase transitions, of classical Heisenberg spins with ex¬ 
change and dipolar interactions on the frustrated kagome 
lattice. 

Our first central result is a determination of the ground 
state for classical Heisenberg dipoles. This is a ferrimag¬ 
netic three-sublattice one. Note that dipolar interactions 
for Heisenberg spins lead to ground states in two dimen¬ 
sional systems effectively confined in the plane of the lat¬ 
tice as a result of e xtens ive energy cost of any finite out 
of plane comp onenfl^R^. In our studies we indeed ob¬ 
serve only in-plane spin-states as the ground states of the 
model. Therefore, the ground states we find also apply 
to classical XY spins in the plane of the lattice. 

Next, we observe that switching on a weak dipolar in¬ 
teraction lifts the extensive ground-state degeneracy of 
the nearest-neighbour model which exists here as it does 
in many other frustrated lattices, e. g. th e Archimedean 
pyrochlore lattice in three dimensions^SBZI^ jg both cases, 
the elementary simplices - triangles for kagome, tetrahe- 
dra for pyrochlore - have vanishing total dipole moment 
in the nearest-neighbour ground state; upon adding dipo¬ 
lar interactions, they enter a state where the quadrupole 
moment of each simplex also vanishe^^. However, for 
stronger values of the dipolar interaction, the suppres¬ 
sion of the leading multipole moment no longer seems 
to be favourable. The general principles governing the 
low-energy states on individual cluster^, and how they 
combine to form a large lattice, is an intriguing topic for 
future studies. 

The concomitant line of phase transitions into the fer¬ 
rimagnetic state at dominant D appear to have expo¬ 
nents v and j3 change monotonically with the ratios (3 jv 
and 7 / 1 Z constant. This is known in the context of the 
’’weak universality” hypothesis and often appears in sys¬ 
tems with n-fold anisotropy where exponents appear to 
‘drift’ from KT values to those of discrete continuous 
transitions. The presence of an enigmatic slice of the 
phase diagram where our methods fail to produce a re¬ 
liable answer further focuses attention on the possibility 
of the appearance of incommensurate states for delicately 
balanced exchange and dipolar interactions. 

Moreover, it seems rather remarkable that the flat 
band of zero-energy excitations simply moves up in en¬ 
ergy without acquiring almost any dispersion. We do 
note that this phenomenon is not so uncommon, after 
all, with a range of different perturbations capable of 
producing a similar phenomenon, a case in point being 
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magnetoelastic interaction^^. Also, a recent preprint!^ 
noted the same phenomenon for a dipolar magnet on the 
Gadolinium Gallium Garnet (GGG) lattice, which has 
historically played an immensely important role in the ex¬ 
perimental study of frustrated magnetic materials. This 
may very well be one of the best experimental handles on 
dipolar interactions, leading to an almost k-in dependent 
resonance in inelastic neutron scatterinj2lllll at a non¬ 
zero energy scaling quite sensitively with the size of the 
dipolar interaction ^ \/~D. 

The prospect for experimental work in this field is 
probably better now than it has been for a very long 
time. In an large number of systems the role o f dipolar 
interactions is important or even dominaniPISll There is 
significant progress in fabrication of dipolar nano arrays 
with a complex frustrated lattice geometrjl^lllil as well 
as recent progr ess on building a dipolar systems in opti¬ 
cal lattice^^*^. In addition recent prog ress in fabricating 
thin films of frustrated materials^^^ suggests a possi¬ 
ble route for realization of dipolar films with a kagome 
geometry. Here the possible candidates for a film real- 
izatio n cou ld be fee kagome materials RhMn 3 , PtMn 3 , 
IrMnqEHHl whe re the latter one is commonly used in thin 
film technologj^^^Ml^ RhMn 3 and PtMn 3 have the fee 
crystal structur^^lHsll where magnetic Mn ions reside on 
the cube faces and the nonmagnetic (Ir) ions site at the 
cube corners. The magnetic ions can thus be viewed 
as being on ABC stacked (111) kagome planes, where 
each site has eight NNs (four in-plane, two to the plane 
above, and two to the plane below). The (111) plane 
is perpendicular to the film plane in thin-film applica¬ 
tions and thus one deals with a thin stack of L kagome 
layers. Interestingly the bulk of IrMn 3 exhibits a long- 
range magnetic order below Tat « 960which is the 
3D manifestion of the 120° q = 0 spin structur^^ one 
of the structures found in our studies to be stabilized by 
weak dipolar interactions. Similar magnetic order is also 
found in RhMn 3 and PtMn 3 . 

We hope that our work will provide motivation for de¬ 
tailed characterisation of nature and collective behaviour 
of some of these experimental systems. 
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along its local z-axis so that it is related to the spin 
operator defined in the crystallographic frame via the 
rotation: 


S,(R'=) = K-/S,(R'=) (Al) 


where Ki is the corresponding rotation matrix. In the 
local frame the Hamiltonian reads 

^ E E E (Ry')5“(R'=)5^(R') (A2) 

i,j a,13 k,l 

where the interaction matrix components transform as 
Ji,(R^)=K:ij,,(R“)K,. (A3) 


Fourier transforming spin operators and interaction ma¬ 
trix gives 


^“(R'^) = ^ E (k) exp [^k ■ (R'= -f r,)] (A4) 


J)y(k)=E^(R“)exp 


kl 


—zk-R! 


kl 


(A5) 


where N is the number of underlying Bravais lattice 
points. Thus, the Hamiltonian in reciprocal space is 


^ = -^EEE^f(k)j(^'W‘^^(-k)- (A6) 

i,j a,13 k 


The linearized Holstein-Primakoff transformation then 
gives 




'S'*"'(k) = V 2 [G(k) -kCi(-k) 
^ r 


(k) = *Y ^ [c|(k) - c,(-k) 
Sf (k) = vV SSk^o exp [-ik • r^] 


(A7) 


1 

Vn 


Ecl(k')G(k'-k), 


with boson operators 
only terms up to second order, we obtain 


Ci(k), cj(k') = 4,k'- Keeping 
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Appendix A: Linear spin-wave theory 

Our spin wave analysis in the non-collinear magnetic 
systems starts with a rotation from the global z-direction 
to the local frame for each moment. Let Si(R^) point 


H = -F 77(1) ^(2) 


(A8) 


where 


i/(o)=_lA52^J)f(0) 

77(1) ^ [f.,( 0)4(0) +F*.(0)cl(0)' 

= -^^EE [A.,(k)cl(k)c,(k) +B,,(k)c|(k)cJ.(-k) 

k 

+B*-(k)c,(-k)cj(k) -h A* (k)c,(-k)cj(-k)j (A9) 
























and 


F.,(o) = j,“(o) + *j;f(o) 

A, ,(k) = i {j;r (k) + jyy(k) - ^ [j^y^k) - j,f (k)]} 
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B. 7k) = I {7.r(k) - (k) + * [j;7(k) + J,f (k)]}. 


where (k) and (k) are new boson operators and all 
the eigenenergies ei(k) are real. 

The specific heat is 

[ei(k)nB(ei(k))]^ exp [/3ei(k)] (A12) 

k i 

where nB(ei(k)) = (ei(k) — 1)“^ is a Bose factor. The 
sublattice magnetization M{T) is obtained by taking into 
account the role of quantum and thermal fluctuations: 


The equilibrium condition that on every site the effective 
magnetic field be parallel to the spin direction implies the 
absence of linear terms. This is satisfied if Fij(O) = 0. 
If the spin ground state is stable after the canonical trans¬ 
formation the Hamiltonian can be written in diagonal 
form 


Ff = Ff^°)+EE^*(k) (All) 

k i 


+EE«*(^) 


a|(k)oj(k) -I- a|(-k)aj(-k) 


M(T) = ^-A5-^EE ’^B(e*(k)) (A13) 

k i 

where 

^S=((vSE1Q’Q1..-i) (AM) 

is the zero-temperature reduction of classical spin polar¬ 
ization and Q is the matrix diagonalizing the spin-wave 
Hamiltonian. 
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